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Abstract 

The diffusive non-Markovian motion over a single-well potential barrier in the presence of a 
weak sinusoidal time-modulation is studied. We found non-monotonic dependence of the mean 
escape time from the barrier on a frequency of the periodic modulation that is character to the 
stochastic resonance phenomenon. The resonant acceleration of diffusion over the barrier occurs at 
the frequency inversely proportional to the mean first-passage time for the motion in the absence 
of the time-modulation. 



I. INTRODUCTION 



The response of complex nonlinear systems on periodic time input may have the features 
that are absent for linear systems. The very famous example of such features is the stochastic 
resonance phenomenon l|, when the response of the nonlinear system on the harmonic 
perturbation is resonantly activated under some optimal level of a noise. The resonant 
activation of the system occurs when the frequency of the modulation is near the Kramers' 
escape rate of the transitions from one potential well to another. The stochastic resonance 



phenomenon has 



3een found and studied in many physical systems like a ring laser 
3], optical bistable systems {4] and others. The prototype of the stochastic 



magnetic systems 

resonance studies is a model of overdamped motion between potential wells of the bistable 
system. The frequency of the transitions between wells is given by the famous Kramers' rate 
and the stochastic resonance is achived when a frequency of an external periodic modulation 
is of the order of the Kramers' rate. 

In the present paper, we would like to study a diffusion over a single-well potential barrier 
in the presence of periodic time modulation. The diffusion over the barrier is generated by 
a colored noise whose statistical properties are related to the retarded friction term. In 
addition to the works , we are going to study how the correlation time of the colored 

noise influence the first passage time distribution and escape rate over the barrier. 

The paper is organized as follows. In Sect. HIl we set in basic Langevin equation of motion 
for diffusive motion over a potential barrier in the presence of sinusoidal time modulation. In 
Sect, nil At it is considered the unperturbed path of the model system. The time modulated 
diffusion is discussed in Sect. IIII B[ Finally, the main conclusions of the paper are given in 
Sect. Ml 



II. DIFFUSIVE MOTION OVER A POTENTIAL BARRIER 

We start from a quite general Langevin formulation of the problem of diffusive overcoming 
of a potential barrier in the presence of a harmonic perturbation Vex_t{t) = a ■ sin(a;t): 

Mqit) = - K{t - t')q{t')dt' + ^(t) + a ■ sin(a;t), (1) 
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where q is the dimensionless coordinate, M is the mass, n{t — t') is the memory kernel 
and C,{t) is the random force. The potential energy i^pot is schematically shown in Fig. 1 
and presents a single- well barrier formed by a smoothing joining at q = q* of the potential 
minimum oscillator with the inverted oscillator, 

^pot = ^Muj\{q - qA^, q < q* , 
= ^pot,B - ^Muliq - qB)\ q > q* . (2) 

A noise term ^(t) in Eq. ([1]) is assumed to be Gaussian distributed with zero mean and 
correlation function related to a memory kernel K{t — t') of a retarded friction force: 

mm) = TK{t-t'). (3) 

Below we shall assume that the memory kernel is given by 

K{t - t') = Koexp (^-^-^—^^ , (4) 

where r is a correlation time. 



III. NUMERICAL CALCULATIONS 



In the numerical calculations, we have measured all quantities quantities of the dimension 
of energy in units of the temperature of the system Eq = T, quantities of the dimension 
of time in units of to = y^M/T and quantities of the dimension of frequency in units of 
Uq = ^T/M. For the system's parameters g^, q*, qB,ujA, ujb, -Spot,B and kq, we adopted 
the values: 

qA = 1, q* = 1.2, qs = 1.6, 
oja = 6.75, ub = 9.59, Epot.B = 5.15, kq = 1920. (5) 

which are widely used under model diffusion-like studies of fission of highly excited atomic 
nuclei, see Ref. jsl]. 
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A. Unperturbed diffusion over the barrier 



At the beginning, we investigated the non-Markovian diffusive dynamics for the infinitely 
slow {u = 0) modulation and calculated a first-passage time distribution. For that, the 
Langevin equation ([T]) was solved numerically by generating a bunch of the trajectories, all 
starting at the potential well (point A in Fig. 1) and having the initial velocities distributed 
according to the Maxwell-Boltzman distribution. First, we would like to study the diffusive 
dynamics in terms of first passage time distribution. 

In Fig. 2, the mean first-passage time is presented as a function of the correlation time 

T. 

An increase of the mean first-passage time tmfpt with the relaxation time r means that 
memory effects in the Langevin dynamics ([1]) hinders the diffusion over the barrier. A non- 
monotonic growth of the mean first-passage time is caused by transition (occurred at the 
T ^ 0.007 ) from nearly Markovian regime of the diffusion ([1]) to the regime where the 
memory effects are quite important. As far as the memory effects become stronger and 
stronger (r — t- oo), the value of tmfpt reaches a finite limit. Here, the hindrance of the escape 
over the barrier is caused by the renormalization of the ordinary conservative force in Eq. (JT] 
) that obtains an additional contribution from the time-retarded force 

- Ko ^ exp (^-^-^—^"^ q{t')dt' -> -kq [q{t) -Qa], t ^ oo. (6) 

In the opposite limit of quite small values of the correlation time r, the hindrance is 
exclusively due to an usual friction: 

- Ko^ exp (^-^-^—^^ q{t')dt' -K^Tqit), r ^ 0. (7) 

In intermediate case, the hindrance is defined both by the friction and the additional elastic 
force, 

- Ko f exp f ) q{t') dt' = -^{t, r) q{t) - C{t, r) g(t), (8) 



/o V , 

where the effective friction coefficient ^{t, r) may be quite well approximated by 
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The stiffness parameter C{t,T) in Eq. ([8]) grows from to kq with a growing of r. 

Secondly, we measure the diffusion dynamics through an escape rate, R{t), char- 

acteristics. The escape rate over the barrier is defined as follows 



m)- ^ ^^^^^ no) 

where P{t) is the survival probability, i.e., probability of finding the system on the left from 
the top of the barrier up to time t: 

Pit) = N{t)/No. (11) 

Here, N{t) is a number of the trajectories being not reaching the top of the barrier up to 
time t and A^o is a total number of the trajectories involved in the calculations. In Fig. 3, 
it is plotted the typical time behavior of the escape rate R{t) for quite small r = 0.005 and 
fairly large r = 0.026 values of the correlation time. 

It is seen from Fig. 3 that initially the escape is affected by transient effects, when 
the survival probability P{t) deviates strongly from the exponential form. With time, the 
escape process becomes more and more stationary giving rise to the corresponding saturation 
of the rate R{t) (JTOl) establishing of a quasistationary probability flow over the barrier. 
Qualitatively, one can describe typical time evolution of the escape rate as 

i?(t) = i?o(l -e"*/*"-^"). (12) 

In both cases a duration of the transient period, ttran, is almost the same (ttran ~ 50) for 
quite weak and fairly large memory effects in the diffusion process. However, a saturation 
value, Rq, of the escape rate is significantly different. It is because of the memory effect for 
the large values of the correlation time r. 

In Fig. 4, we showed how the value Rq f|T2!) depends on the size of the memory effects in 
the diffusive dynamics ([T]). Dotted line in Fig. 4 represents the famous Kramers' result for 
the escape rate value in a quasistationary regime jsj 



1 + 



7(r) 



2MuJi 



> exp f-l ) , (13) 



2Mur ' \ T 



where the r-dependent friction coefficient 7(r) is given by Eq. 
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We see that the memory effects significantly suppress the value of the escape rate in the 
saturation regime of probability flow over the potential barrier. Initially (i. e., at relatively 
small values of the correlation time r) the suppression is mainly caused by the growing role 
of the usual friction in the non-Markovian diffusion ([T])- see also Eqs. ([H]) and As 
is followed from the Fig. 4, in this case the escape rate at saturation Rq (IT^ may be quite 
well approximated by the Kramers' formula (fT3l) . On the other hand, at relatively large 
correlation times r, the effect of the friction on the diffusion over the barrier is negligibly 
weak and the escape rate's suppression appears exclusively due to the additional conservative 
force, see Eq. ([HD. As a result, the stationary value of the escape rate deviates substantially 
from the Kramers' escape rate fll3p at the fairly strong memory effects in the diffusive motion 
across the barrier. 

B. Diffusion over the barrier in the presence of a periodic modulation 

Now we will study the diffusion over the barrier in the presence of the external 

harmonic force. We will assume that the amplitude a of the force a ■ sin(co't) in Eq. ([1]) is 
so small {a = 0.05) that still the reaching the top of the barrier is caused exclusively by 
diffusive nature of the process. 

In Fig. 5, we calculated the typical dependencies of the mean first-passage time tmfpt on 
the frequency u of the external harmonic force. The calculations were performed for the 
weak, T = 0.005, (lower curve in Fig. 5) and strong, r = 0.026, (upper curve in Fig. 5) 
memory effect on the non-Markovian diffusive motion over the barrier. 

In both cases the mean first-passage time tmfpt non-monotonically depends on the fre- 
quency of the perturbation that is character for the stochastic resonance phenomenon ob- 
served in a number of different physical systems. From Fig. 5 one can conclude that diffusion 
over the potential barrier in the presence of the harmonic time perturbation is maximally 
accelerated at some definite so to say resonant frequency u^es of the perturbation. 



see also Figs. 2 and 4. In fact, the quantity tmfpt(i^ = 0) presents the characteristic time 
scale for the diffusion dynamics ([I]). In the case of adiabatically slow time variations of the 
harmonic force, c<;tmfpt(i^ = 0) ^ 1 and t < tmfpt, one can approximately use a-sm{cot) ~ a-ut 



1.5 



(14) 



res 
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and the diffusion over the barrier is shghtly accelerated. As a result of that, the mean first- 
passage time tmfpt(i^) is smaller than the corresponding unperturbed value tmfpt(i^ = 0). The 
same feature is also observed at the fairly large modulation's frequencies. Thus, in the case 
of ut^{pt{u = 0) ^ 1, the harmonic perturbation a ■ sin(a;t) may be treated as a random 
noise term with the zero mean value and variance a^. Such a new stochastic term will lead 
to additional acceleration of the diffusion over the barrier. 

The existence of the resonant regime f|T^ in the periodically modulated diffusion process 
(II]) is even more clear visible in the escape rate characteristics of the process. We have 
plotted in Fig. 6 the time evolution of the escape rate f lTU]) found for the resonant frequency 
i^res dHl) (curve 1 in Fig. 6), the quite smaller frequency u = Wres/lO (curve 2 in the Fig. 6) 
and the quite larger frequency u = 10 Wres (curve 3 in Fig. 6) of the modulation. 

Again, at very slow (w = Wros/lO) and fast (w = lOwres) perturbations the escape rate R{t) 
looks very similar to the corresponding unperturbed value R{t, u = 0), compare Figs. 3 and 6. 
In other words, the initial transient period in the time evolution of the escape rate is followed 
by the stationary regime, when the escape subsequently saturates with time. Contrary to 
that, the escape rate shows complex time behavior as far as the periodic modulation occurs 
at the resonant frequency {u = Wrcs), see curve 1 in Fig. 6. This resonant regime of the 
modulated diffusion is essentially nonstationary when the system remains excited during 
quite long time. We checked such a feature for larger number of trajectories and longer time 
intervals used to calculate the escape rate characteristics ( ITOl) . 

IV. CONCLUSIONS 

In the present study, we have investigated how model non-Markovian diffusion over the 
single-well parabolic barrier is affected by the external periodic time modulation. We have 
calculated both the mean first-passage time tmfpt and the escape rate R(t) over the barrier. 
These two quantities have been found to be sensitive to the relative size of memory effects 
in the diffusive dynamics measured by the correlation time r. Thus, we have 

demonstrated that the memory effects hinder the escape over the barrier, see Figs. 2 and 
4. In contrast to the motion in the presence of usual friction force, the hindrance of the 
escape occurs due to the Markovian friction and additional conservative components ([H]) of 
the retarded time force in Eq. ([T]). Having calculated the mean first-passage time tmfpt for 
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different values of the frequency u of the modulation, we have found that the sinusoidal 
perturbation accelerates the diffusion over the barrier, see Fig. 5. The maximal (resonant) 
acceleration is achieved at the u = Uj-es, where u^es is inversely proportional to the mean 
first-passage time in the absence of the modulation f[T^ . We have seen that a value of the 
resonant activation over the barrier tmfpt(i^ = i^rcs)/^mfpt(i^ = 0) remains practically the same 
for the quite weak as well as for the fairly strong memory effects in the diffusive dynamics. It 
has been observed that the diffusive dynamics in the resonant activation regime f lT^ 

has the peculiarity, reflecting in the complex time behavior of the escape rate R{t), see Fig. 6. 
Importantly, that the absence of the escape rate's saturation with time implies essentially 
transient character of the events of the first passing the top of the potential barrier. 
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Figure captions. 



Fig. 1. Landscape of potential energy Epot{q) of Eq. ([2]). 

Fig. 2. The mean first-passage time tmfpt of the non-Markovian diffusion process 
vs the correlation time r. 

Fig. 3. The time dependence of the escape rate R{t) (fTO!) of the non-Markovian diffusion 
process ([1])- (jl]) calculated for the cases of quite small r = 0.005 and fairly large r = 0.026 
values of the correlation time. 

Fig. 4. The saturation value Rq of the escape rate f|T0|) - f|T2l) vs the strength r of the mem- 
ory effects in the non-Markovian diffusion process The dotted line represents the 
Kramers' result f[T^ for the escape rate calculated with the r-dependent friction coefficient 
ofEq. ©. 

Fig. 5. The mean first-passage time tmfpt of the non-Markovian diffusion process ([I])-(jl]) 
is given as a function of the frequency u of the harmonic time perturbation at two values of 
the correlation time r = 0.005 (lower curve) and r = 0.026 (upper curve). 

Fig. 6. The time dependence of the escape rate R ffTOl) for the periodically modulated 
diffusion over the barrier ([l])-(j4]) are shown at the resonant frequency cjres (fT4l) (curve 1), 
fairly smaller frequency u = Wrcs/lO (curve 2) and quite bigger frequency u = lOwrcs (curve 
3) of the periodic modulation. The dependencies were calculated for the memory time 
r = 0.01. 
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